clear
close all
clc
delete *.asv

cd ..\
cd ..\
cd Stored_File

load Comparative_kappa
load Comparative_epsilon
load Comparative_betadelta

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%% FIGURE C1 %%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Ross measure for risk
subplot(3,3,1)
set(0,'DefaultAxesTitleFontWeight','normal')
plot(kappa_vector,(Ross_kappa),'k','LineWidth',0.5)
hold on
set(gca,'FontSize',05);
ylabel('$${\Upsilon}$$','Interpreter','Latex','FontSize',08)
xlabel('\kappa','FontSize',08)
set(gca, 'Xtick', [0.6 0.625 0.65], 'Ytick', [0.05 0.37], 'box', 'off');
axis([min(kappa_vector),0.65,min(Ross_kappa),max(Ross_kappa)])
box on

subplot(3,3,2)
set(0,'DefaultAxesTitleFontWeight','normal')
plot(sigma_vector,(Ross_sigma),'k','LineWidth',0.5)
hold on
set(gca,'FontSize',05);
xlabel('\epsilon','FontSize',08)
set(gca, 'Xtick', [0.04 0.07 0.10], 'Ytick', [0.08 0.37], 'box', 'off');
axis([min(sigma_vector),0.10,0.08,max(Ross_sigma)])
box on

subplot(3,3,3)
set(0,'DefaultAxesTitleFontWeight','normal')
plot(delta_vector,(Ross_betadelta),'k','LineWidth',0.5)
hold on
set(gca,'FontSize',05);
xlabel('\delta=\beta','FontSize',08)
set(gca, 'Xtick', [0.65 0.82 0.98], 'Ytick', [0.35 1.41], 'box', 'off');
axis([0.65,0.987,min(Ross_betadelta),max(Ross_betadelta)])
box on

% Equivalent Consumption Welfare
subplot(3,3,4)
plot(kappa_vector,1-IC_kappa,'k','LineWidth',0.5)
set(gca,'FontSize',05);
ylabel('$$1-E_{\phi }[\iota (x)] $$','Interpreter','Latex','FontSize',08)
xlabel('\kappa','FontSize',08)
axis([min(kappa_vector),max(kappa_vector),1-max(IC_kappa),1-min(IC_kappa)])
set(gca, 'Xtick', [0.6 0.625 0.65], 'Ytick', [0.05  0.31], 'box', 'off');
box on

subplot(3,3,5)
plot(sigma_vector,1-IC_sigma,'k','LineWidth',0.5)
set(gca,'FontSize',05);
xlabel('\epsilon','FontSize',08)
axis([min(sigma_vector),0.10,0.18,1-min(IC_sigma)])
set(gca, 'Xtick', [0.04 0.07 0.10], 'Ytick', [0.18 0.31], 'box', 'off');
box on

subplot(3,3,6)
plot(delta_vector,1-IC_betadelta,'k','LineWidth',0.5)
set(gca,'FontSize',05);
xlabel('\delta=\beta','FontSize',08)
set(gca, 'Xtick', [0.65 0.82 0.98], 'Ytick', [0.33 0.98], 'box', 'off');
axis([0.65,0.987,min(1-IC_betadelta),max(1-IC_betadelta)])
box on

% Certainty equivalent consumption
subplot(3,3,7)
plot(kappa_vector,CE_kappa,'k','LineWidth',0.5)
set(gca,'FontSize',05);
ylabel('$$E_{\phi }[\theta (\omega )] $$','Interpreter','Latex','FontSize',08)
xlabel('\kappa','FontSize',08)
axis([min(kappa_vector),0.65,min(CE_kappa),max(CE_kappa)])
set(gca, 'Xtick', [0.6 0.625 0.65], 'Ytick', [0.01 0.23], 'box', 'off');
box on

subplot(3,3,8)
plot(sigma_vector,CE_sigma,'k','LineWidth',0.5)
set(gca,'FontSize',05);
xlabel('\epsilon','FontSize',08)
axis([min(sigma_vector),0.10,min(CE_sigma),max(CE_sigma)])
set(gca, 'Xtick', [0.04 0.07 0.10], 'Ytick', [0.02 0.23], 'box', 'off');
box on

subplot(3,3,9)
plot(delta_vector,CE_betadelta,'k','LineWidth',0.5)
set(gca,'FontSize',05);
xlabel('\delta=\beta','FontSize',08)
set(gca, 'Xtick', [0.65 0.82 0.98], 'Ytick', [0.28 2.84], 'box', 'off');
axis([0.65,0.987,min(CE_betadelta),max(CE_betadelta)])
box on
